Shortcuts to adiabaticity: fast-forward approach 
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The "fast-forward" approach by Masuda and Nakamura generates driving potentials to accelerate 
slow quantum adiabatic dynamics. First we present a streamlined version of the formalism that 
produces the main results in a few steps. Then we show the connection between this approach and 
inverse engineering based on Lewis-Riesenfeld invariants. We identify in this manner applications in 
which the engineered potential does not depend on the initial state. Finally we discuss more general 
applications exemplified by wave splitting processes. 
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I. INTRODUCTION 

Motivated by the practical need to accelerate quantum 
adiabatic processes in different contexts (transport iHa] , 
expansions H, Hi , population inversion and control ( 8-fl3l . 
cooling cycles [y, [3, EH > wavefunction splitting flft - 19j |) , 
and by related fundamental questions (about the quan- 
tum limits to the speed of processes, the viability of adia- 
batic computing [20| , or the third principle of thermody- 
namics [14l. l2l| ), a flurry of theoretical and experimental 
activity has been triggered by the proposal of several ap- 
proaches to design "shortcuts to adiabaticity". Among 
other approaches let us mention (i) a transitionless track- 
ing algorithm or "counterdiabatic" approach that adds 
to the original Hamiltonian extra terms to cancel transi- 
tions in the adiabatic or superadiabatic bases [ gjR| ; (ii) 
inverse engineering of the external drivin g [3] , U, |g, l2ll - 
[26j based on Lewis- Riesenfeldt invariants |27I |T which has 
been applied in several expansion experimen ts |25l . |26| ; 
(iii) optimal control (OC) methods [a, 0, EH, Il6l|. some- 
times combined with other methods to enhance their per- 
formance @, H, 0]; (iv) the fast- forward (FF) approach 
advocated by Masuda and Nakamura dHHH; (v) paral- 
lel adiabatic passage [291432} . 

The multiplicity of approaches is quite useful because 
they may complement each other: either in the same ap- 
plication, as demonstrated e.g. with OC and invariant- 
based methods, or because of their different domains. 
Clarifying the features, overlaps, and relations among 
these approaches is important to apply the ones which 
are best suited for specific systems and objectives [33| . 
or to develop new ones. In this paper we shall establish 
in particular the connection between the fast-forward and 
the invariant-based methods. 

Based on some earlier results [34| , the fast- forward for- 
malism for adiabatic dynamics and several application 
examples were worked out in [l9|, [28| by Masuda and 
Nakamura for the Gross-Pitaevskii (CP) or the corre- 
sponding Schrddinger equations. The objective of the 
method is to accelerate a "standard" system subjected 
to a slow variation of external parameters. The time is 
then rescaled by a "magnification factor" , and an ansatz 
wave function is defined by the standard function mul- 



tiplied by a phase factor that, in general, depends on 
position and time. Inserting the ansatz into the dynam- 
ical equation provides in principle the form of the neces- 
sary fast-forward driving potential and the equation to 
be satisfied by the phase. This procedure leads however 
to divergent terms when the reference standard process is 
infinitely slow. The solution found in |19j to this problem 
was to "regularize" the Hamiltonian and standard state 
using a new ansatz. 

As a consequence of the different steps and functions 
introduced the resulting procedure is somewhat involved, 
which hinders a broader application. In Sec. [II] we pro- 
vide a streamlined construction of local and real fast- 
forward potentials, and Sec. Mil delves into a more de- 
tailed connection between this streamlined version and 
the original formulation of the FF formalism. Section [TV] 
relates the fast-forward approach to the inverse method 
for dynamical invariants which are quadratic in momen- 
tum. Section [V] discusses applications beyond this do- 
main, in particular wavefunction splitting, which is an 
important operation for matter wave interferometry [fH - 
El, [35|. Finally Sec. IVII discusses the results and open 
questions. 



II. A SIMPLE INVERSE METHOD 

Our starting point is the 3D time-dependent CP equa- 
tion 



(i) 



where the Hamiltonian H is the sum of the kinetic energy 
T, the external potential V(t), and the mean field poten- 
tial G(t). We are assuming an external local potential, 
where "local" means here (x.\V(t)\x.') = V(x, t)<5(x — x'). 
Then, by solving Eq. ([I]) in coordinate space, V(x,t) 
may be written as 



V(x,i) 



(2) 



2 



with (x\ip(t)) = ip(x, t). The kinetic and mean field terms 
in the coordinate representation have the usual forms 



III. CONNECTION WITH THE 
FAST-FORWARD APPROACH 



mm) 



2m 



V 2 V(x,*), 



<x|G(i) !</,(*)) = g\^(x,t)\ 2 i>(x,t), 

g being the coupling constant of the Bose-Einstein con- 
densate. By introducing into Eq. ([2]) the ansatz 

(x|</>(*)) = r(x,t)e^ x < 4 \ r(x,t),0(x,t) € K, (3) 
we get 

ft 2 /2iV0-Vr , 



V(x,i) = z/i ^ + 

r 2m 



V 2 r 



(4) 



where the dot means time derivative. The real and imag- 
inary parts are 

Re[V(x,t)] = -h<P+-^(— -(Vcb) 2 ) -gr 2 ,(5) 
2m V r I 



im[y(x,t)] = h- + 

r Im 



r H 2 /2V</> ■ Vr 



(6) 



Our purpose now is to design a local and real poten- 
tial such that an initial eigenstate of the initial Hamilto- 
nian, typically the ground state but it could be otherwise, 
evolves in a time tf into the corresponding eigenstate of 
the final Hamiltonian (a different goal will be discussed 
in the final Section). We assume that the full Hamilto- 
nian and the corresponding eigenstates are known at the 
boundary times. 

By construction the potential of Eq. (j4|) is local. If we 
impose Im[V(x, t)] = 0, i.e. 



- + — — + V 2 < 

r Im \ r 



= 0, 



(7) 



then we get from Eq. ([5]) a local and real potential. 

In the inversion protocol we design r(x, t) first, then 
solve for <p in Eq. J7j), and finally get the potential V 
from Eq. ([5]). If, at the boundary times, r = is im- 
posed, Eq. (J?]) has solutions <f>(x, t) fulfilling that </)(x, t) 
is independent of x at t = and t — tf. Using this in 
Eq. (O at t = 0, and multiplying by e l ^°\ we get 



■ — V 2 +U(x,0)+ ff ^(x,0)| 5 
Im 



tp(x,0) = -ft0(O)^(x,O). 

The initial state is thus an eigenstate of the stationary 
GP equation at t = 0, and —h<ft(0) = E(0) is the energy 
of the eigenstate ip(x,0). Note that the above solution 
of <f> (with r — at boundary times) admits the addition 
of an arbitrary function that depends only on time and 
modifies the zero of energy. A similar result is found at 



We shall now reformulate the above results to conect 
them with the FF approach in pjj [34[ ■ The notation is 
made close but not necessarily in full agreement with [Til . 
l34j | . Let us define an external parameter that depends on 
time, or on some scaled time function, according to 



R(A(t)) = eA(i) =: K(t). 



(9) 



Here 7Z(t) and R(A) are in general different functions of 
their arguments, e is a small positive constant and the 
scaling function A(t) is given in terms of a magnification 
factor a, 



A(t) 



dt'a(t') 



(10) 



a(t) is positive for < t < tf and zero at the boundaries 
t = and t = tf. Note that 1Z = ea and It = ea. We 
rewrite the modulus and the phase in Eq. ([3]) as 



r(x,t) = r(x,K(t)), 



dt'£(TZ(t')) 



(11) 

ea(t)e(x,1Z(t)), (12) 



where again we have distinguished the functions accord- 
ing to their different arguments, in particular E(t) = 
£(lZ(t)). If we also demand a — at the bound- 
aries to fulfill r = 0, then 0(x,O) = -8(11(0)) /h and 
(£(x,*/) = -£(H(t f ))/%. Substituting Eq. (jlS) in Eq. 
(J7J, 9 has to satisfy 

= f(x,K(t))\7 2 9(x,K(t)) + 2Vf(x,ft(i)) ■ V9(x,TZ(t)) 



2m df , _ . . . 



(13) 



and from Eq. ([5]), the "fast- forward" potential is given 
by 

V(x,t) = V (x,TZ(t)) - hea(t)6(x,n(t)) 
d6 



-he 2 a 2 (t)^(x,TZ(t)) 
-^a 2 (t)[Vd(x,K(t))} 



(14) 



where the "standard potential" Vq = Vo(x,1Z) is defined 
by the stationary GP equation 



V 2 + Vb(x,ft)+ 3 f 2 (x,7e) 

Im 



r(x,TZ) = S(R.)r(x,K). 

At the boundary times, but in general only there, 
V(x,t)=V Q (x,TZ(t)). 

Equation (|14l) for the driving potential coincides with 
Eq. (2.28) in 19] for real eigenfunctions, whereas Eq. 
(IT31 for the phase function 9 corresponds to Eq. (2.18) 
in [l|. 
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The present formal framework may be used in the fol- 
lowing way: (i) Starting from a given standard potential 
Vb(x,K), r(x,7e) and £{K) would follow from Eq. fig]) . 
Alternatively, it is also possible to impose r(x, 7?.) first 
and then calculate Vo- (ii) An auxiliary function lZ(t) 
is imposed, (hi) 6 has to be determined from Eq. (|13jl. 
(iv) The fast-forward potential can be calculated from 

Eq. mm. 

To arrive at this recipe in [l9l l34j preliminary steps 
are the definition of a standard state, a virtually fast- 
forwarded state, and a regularized state with their cor- 
responding equations. The route followed in Sec. ITTI to 
the driving potential is in comparison quite direct. This 
is so because we made no explicit use of a slow reference 
adiabatic process, although it might be deduced from the 
fast designed dynamics if required. The key simplifica- 
tion is to start with the ansatz in Eq. (J3)) and derive 
the two basic equations for phase and potential from it 
by imposing locality and reality of the driving potential. 
Since the phase <f) that solves Eq. ([7]) depends in general 
on the particular r(x, t), the potential calculated through 
Eqs. (JS]) resp. (fT4"]) gives in principle a state-dependent 
potential. However, in some special circumstances, as we 
shall see below, the fast-forward potential becomes state 
independent. 



and a = a(t) must satisfy the auxiliary equations 

/ .2 



F(i) 



= a 



(17) 
(18) 



with ujq an arbitrary constant. The associated dynamical 
invariant, up to a constant factor, is given by 



I = 



1 

2m 
1 



|p(p — ma) — mp(x — a)| z 



-mu>t\<T\ z + U(<T), 



(19) 



with p = —ihV. It satisfies dl/dt — dl(t)/dt — 
i [I(t), H(t)] = 0, so its expectation values are constant 
for any wave function ip(t) that evolves with H. 

For the potential in Eq. (fl~6|) . the general solution of 
the time-dependent Schrodinger equation, Eq. ([T]), can 
be expanded as a linear combination with constant coef- 
ficients c n and orthonormal eigenvectors ipn of I [13] , 

tf(x,t) = ]T C „e M »^(x,i), (20) 

n 

/Vn(x,i) - X n ^ n {x,t), (21) 



IV. CONNECTION WITH INVARIANT'S 
BASED INVERSE ENGINEERING APPROACH 



In this section we shall relate the previous results for 
the linear (g = 0) Schrodinger equation to the engineer- 
ing approach based on quadratic-in-momcntum invari- 
ants. The non-linear GP equation could also be treated, 
as in [3, [22j], but it does not allow in general for the 
state-independent potential forms that we shall describe 
for 5 = 0. 



A. Lewis-Leach potentials 

In a direct (rather than inverse) approach, the po- 
tential V(x, t) is considered to be known, and the wave 
function at any time t can be deduced from the Lewis- 
Riesenfeld theory of invariants [27]. Suppose that the 
potential V(x, t) has the structure of the most gen- 
eral "Lewis-Leach" potential that admits a quadratic-in- 
momentum invariant 13611 . 



V(x,t) = -F(t)-x+lmLu 2 {t)\x\ 2 + \u(cr) + h(t), (16) 
2 p £ 



where u)(t), F(t) and h(t) are arbitrary functions of time 
and U (<t) is an arbitrary function of its argument er = 
<x(t) = (x— a)/ p. The time dependent functions p = p(t) 



where A„ are the time independent eigenvalues of /. The 
phases atn satisfy h^f- = (t/j n \ih§- t - H\^ n ), 0,3^3 



dt' 



A, 



m[\ap — ap\ 



-^o 2 I«I7p 2 ] 



Performing now the unitary transformation 0, 0] 

1 



-h 
(22) 



^,(x,t) = efl^W /^+(^-ap)-/p]_ Xii((T)i (23) 



the state ipn is easily obtained from the solution Xn(°") 
(normalized in cr-space) of the auxiliary stationary 
Schrodinger equation 



h 2 2 
— 

2m 



-U(cr) 



X n {tr) = KXn{(T). (24) 



In the direct approach we assume that U(cr), u>(t) and 
F(t) are known. Solving Eqs. (JT7J! and HU) we get p(t) 
and a(t) from them. Thus Eq. (|24p can be solved to get 
A„ and Xn(cr). Finally combining Eqs. (|23|) and ([24]) . 
the mode e la ' l ipn can be calculated at any time. 



B. Inverse engineering approach 

In the inverse approach based on quadratic-in- 
momentum invariants, the Hamiltonian is assumed to 
have the form given in Eq. (|16[) . at all times and in 
particular at initial and final instants. As U is given the 
stationary Eq. (j2"4"]l may be solved. Then the functions 



4 



p and a are designed so that [H(t),I(t)] = for t = 
and t = tf. Thus the Hamiltonian and the invariant 
have common eigenvectors at these boundary times 0- 
0, HI, H3] - Typically the initial state tp(0) is the ground 
state of -ff(O) which is also an eigenstate of 1(0), and this 
state evolves according to Eq. ([20)) . as an eigenvector of 
the invariant. 

To relate the above to the simple inverse method of 
Sec. |TT]we consider the single mode wave function ip n of 
Eq. ((231) and identify 



r„(x,i) = X n(<r)/p 3/2 (t)- 



(25) 



The subscript n underlines the dependence with the nth 
mode considered. Note that p and a are chosen at this 
point. We may get the phase 4> n from Eq. ([7]). It can be 
checked by direct substitution that 



-[p\x.\ 2 /2p+(ap-ap)-xL/p] 



(26) 



where 



Ht) 



dt' 



m[\ap — ap\" 



2p 2 



(27) 

is a solution of this equation. Once r n and the phase <p n 
are known, Eq. §5§ gives the potential V n (x,t). A dif- 
ferent arbitrary function of time J-(t) in (j) n (x.,t) would 
produce a shift of the zero of energy in the resulting po- 
tential V n (x,t). We get, using Eq. ([24"]). 



a ; I x 



m ( ujq + p 3 p 



;U(<t) 



Taking Eqs. ([T7]) and (fT8"|) into account, this potential 
agrees with the potential in Eq. (fT6)) . It is by construc- 
tion local and real. Moreover it is independent of the 
nth state considered so that linear combinations of the 
modes at t = end up at tf unexcited, preserving the 
initial populations. 

Summarizing, the inverse engineering approach and 
the fast-forward approach are connected via the simple 
inversion method. Clearly, the external parameter 7Z in- 
troduced in Sec. Ill must be related to a and p, as 
illustrated in the following section. 



C. Example: harmonic expansion 

Now we discuss an example of a 3D harmonic ex- 
pansion produced with the inverse engineering approach 
based on invariants and with the fast forward technique 
to illustrate the links between the two methods. 



Invariants based approach: Suppose that the expansion 
is governed by the Hamiltonian 



2rn 2 



(28) 



where ui(t) is unknown, but at the boundary times w(0) = 
wo and w(i/) = u)f. This potential is a particular case of 
Eq. (HU) with F(t) = U{cr) = h(t) = 0. Equation QJJ) 
is trivially fulfilled if a(t) = a(t) = and consequently 
ct(t) = 0. p(t) has to satisfy the Ermakov equation, 
Eq. (|17p. The inverse engineering consists on imposing 
conditions on p and its derivatives, 

p(0) = 1, p(0) = 0, p(0) = 0, 
p(t f ) = 7 , p(t f ) = 0, p(t f ) = 0, (29) 

where 7 = (ujo/uif) 1 / 2 , to guarantee the commutation 
between H(t) and I(t) at t — and tf, and then getting 
u)(t) from Eq. ([Tf]). 

Fast forward approach: The starting point for the fast 
forward approach could be the nth eigenstate of a har- 
monic trap (see also [lpj) with angular frequency 1Z = u, 



X„(x,7e) = 



tt 3 /4 y/2 n *+ n y+ n *n x \n y \n z \ 



H nx (px)H ny ((3y) 
(30) 



where f3 — \J mlZ/h. This state plays the role of r. The 
corresponding potential Vq is clearly 



V (x,K) = yi? 2 x 2 , 



and 



£ n = huji n x + n y + n z + 



(31) 



(32) 



The first step is to solve Eq. (fT3"]) and we get as a solution 



9(x,K) = -■ 



&H1Z 



(33) 



Connection: The connection between the auxiliary 
variable 1Z in the fast forward approach and the aux- 
iliary variable p(t) in the inverse engineering approach is 
in this example explicitely given by 



n{t) 



^0 
p(t) 2 



(34) 



see Eqs. (|2"5|) and The boundary conditions for p(t) 
in Eq. (|2^|) become 

n(o) = uj q , n(o) = 0, ii(o) = o, 
n{t f ) = Uf , ii(t f ) = 0, ii{t f ) = 0. (35) 

It also follows that ea(t) = K(t) = -2u p(t) / p(t) 3 . The 
auxiliary functions p(t) resp. TZ(t) can be chosen in some 
way fulfilling the boundary conditions. 
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The corresponding potential in the inverse engineering 
formalism is constructed by first solving the Ermakov 
equation to get ui 2 (t). Then one has 

V-l^WW-lm^-^W, (36) 

whereas the fast-forward potential is given, according to 
Eq. flU]), by 

h 2 2 n ,2 m 2 A 
2m fc U 4/i 2 7?, 2 y 

_ m|x| 2 / _ p(t) \ 

- 2 ^ p(t)* p^j' 

which agrees with Eq. (|36|) . 




x (urn) 



FIG. 1: (Color online) Phase 4>(x,t) calculated from Eq. 
Parameters: m = 1.44 ■ 1CT 25 kg, a = 3/im, oj/27t = 125 Hz 
and i/ = 80 ms. 



V. BEYOND LEWIS-LEACH POTENTIALS: 
WAVEFUNCTION SPLITTING PROCESSES 

The transitionless condition for the inverse engineering 
method based on invariants relies on the commutativity 
[H(t),I(t)] — at times t = and tf, which guaran- 
tees common eigenvectors for H and / at these boundary 
times. According to Eq. (|23)l the structure of the density 
of the nth mode of the invariant at initial and final times 
for quadratic-in-p invariants is 



1 



P 3 (0) 



Xn 



a(0) 



P(0) 



1 



P 3 (tf) 



Xn 



'(*/) 



p(V) 



(37) 

The final density is a translation and/or scaling of the 
initial one. This means that for processes in which the 
initial and final eigenstates of the Hamiltonian do not 
behave according to Eq. (|37l) . the commutativity of H 
and I at the boundary times cannot be achieved. This 
restriction is due to the use of quadratic-in-p invariants, 
not to the invariants-based method. Studying and ap- 
plying more general invariants is still an open question. 

As an example in which Eq. (|37|) does not hold for 
the final densities, let us consider the splitting of an ini- 
tial state from a single to a double well potential. For 
simplicity we take the ID linear Schrodinger equation 
governed by the Hamiltonian 



H(t) 



P 

2m 



1 



moj 2 (t)x 2 + rj(t)x 4 



(38) 



For the initial single trap we consider u> 2 (0) = u> 2 and 
77(0) = rjo > 0. The final double well is characterized by a 



repulsive harmonic part with uj (tf) 



id nitf) 



rjf > 0. Comparing terms with Eq. (fT6|) , F = h = a = 
and consequently a = a = 0, U(a) = rj(t)p 2 x i , and the 
particular structure of U(a) sets rj(t) = n/p 6 , where p 
satisfies the Ermakov equation, Eq. (fl~7|) . and k is an 
arbitrary constant. The associated invariant is 



I(t) 



-^—{pp - mpx) 2 + -Too + k- 
2m 2 p z p 



(39) 



Imposing p(0) = 1, p(Q) — 0, p(0) = 0, and identi- 
fying k = 7/0, then H and / commute at t = 0. At 
t f , [H(t f ),I(tf)] = for p(t f ) = ^ /^/) 1/2 , p(t f ) = 
p(tf) = 0, and rjf — —itiquj^/ujq. However, p must be 
a positive real function if initially so, and moreover the 
final potential that we get is complex. 

A way out is to use the simple inverse fast-forward 
method for specific initial and final states without re- 
stricting the potential form, see also [l9[. Consider for 
example the ID splitting of the initial state r(x, 0) — 
e -p x /2 _ m(jJ J ^ j n ^ final form r(x,tf) — 

e -t3\x~af/2 + e -0 2 (x+a) 2 /2^ j n between we apply the 

interpolation 

r{x,t) = z{t)^[l-n(t)]r(x,0)+K{t)r(x,tf)X, (40) 

where lZ(t) is some smooth, monotonously increasing 
function from to 1 and z{t) is a normalization func- 
tion. We also impose that TZ = to ensure r = at 
the boundary times t = and tf. In the numerical ex- 
amples the function lZ(t) is chosen as a polynomial of 
degree 7 to make zero the second and third derivatives 
at the boundaries. Once we have established the form 
of r(x,t), we solve Eq. ([6]), Im[V(x, f)] = 0, to get the 
phase cj> with the initial conditions 4>(0,t) = <p'(0,t) = 
that fix the zero-energy point (the prime means spatial 
derivative). Then the resulting phase is introduced into 
Eq. (JSJ) to get the potential V(x,t). In Fig. [T]the phase 
4>(x, t) is plotted for a non adiabatic process with tf — 80 
ms. The corresponding fast-forward potential is plotted 
in Fig. [5J V(x, t) in Fig. [5] could be realized with high 
resolution time-varying optical potentials "painted" by 
a tightly focused rapidly moving laser beam [371 ] . or by 
means of spatial light modulators [38] . A simpler approx- 
imate approach would involve the combination of three 
Gaussian beams. In principle the time tf can be reduced 
to produce the splitting in a shorter time. For example, 
in Fig. = 10 ms and a more complicated potential 
is needed. 
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FIG. 2: (Color online) Fast-forward potential V(x,t) in units 
of fuo for a final time tf = 80 ms. The rest of parameters are 
the same as in Fig. [T] Note the coexistence of three wells, 
e.g. around 40 ms, before the final two wells are established. 




FIG. 3: (Color online) Fast-forward potential V(x,t) in units 
of fko for a final time tf = 10 ms. The rest of parameters are 
the same as in Fig. Q] 



VI. DISCUSSION 

We have first distilled from the somewhat imposing 
set of equations of the fast-forward (FF) formalism as 



originally presented a streamlined version that may aid 
to apply it more easily. Our second aim has been to 
relate it to other inverse engineering methods. In a pre- 
vious publication, the inverse-engineering method based 
on invariants, was related to the transitionless tracking 
algorithm, and their potential equivalence was demon- 
strated [331 ] . Similarly we have established in this paper 
the connection between the fast-forward method and the 
invariant-based method for quadratic-in-momentum in- 
variants. These relations do not imply the full identity 
of the methods but their overlap and equivalence in a 
common domain. They are still useful heuristically as 
separate approaches since they are formulated in rather 
different terms [lH [33[ . Moreover they facilitate exten- 
sions beyond their common domain, as exemplified by the 
wave-splitting processes discussed in the previous section. 
Further extensions are left for separate analysis: for ex- 
ample the possibility to transfer an excited state into the 
ground state or viceversa, or combining the fast-forward 
approach with optimal control theory (OCT) without in- 
cluding the final fidelity in the cost function as in [ill - fill . 
(This would be possible because the fidelity is guaranteed 
to be one by construction.) It will also be interesting 
for future work to consider complex potentials, either as 
solutions to the shortcut dynamics, as in the quantum 
brachistochrone [39[ , or as an effective description of the 
system dynamics to be accelerated l'2'.j . 
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